International Journal of Power Electronics and Drive System (IJPEDS) 
Vol. 12, No. 1, Mar 2021, pp. 413~420 
ISSN: 2088-8694, DOI: 10.1159 1/1jpeds.v12.11.pp413-420 o 413 


Modeling and simulation of forces applied to the horizontal axis 
wind turbine rotors by the vortex method coupled with the 


method of the blade element 


Ibtissem Barkat!, Abdelouahab Benretem?, Fawaz Massouh’, Issam Meghlaouit, Ahlem Chebel> 
l.2.45Laboratory, Electromechanical Engineering, Badji-Mokhtar-Annaba University, Annaba, Algeria 
ŝLaboratory DynFluid, Arts et Métiers ParisTech, 151 The Boulevard Hospital, Paris, France 


Article Info 
Article history: 


Received Apr 26, 2020 
Revised Jan 14, 2021 
Accepted Feb 4, 2021 


Keywords: 


Forces on the blade 
Method of the blade element 


ABSTRACT 


This article aims to study the forces applied to the rotors of horizontal axis 
wind turbines. The aerodynamics of a turbine are controlled by the flow 
around the rotor, or estimate of air charges on the rotor blades under various 
operating conditions and their relation to the structural dynamics of the rotor 
are critical for design. One of the major challenges in wind turbine 
aerodynamics is to predict the forces on the blade as various methods, 
including blade element moment theory (BEM), the approach that is 
naturally adapted to the simulation of the aerodynamics of wind turbines and 
the dynamic and models (CFD) that describes with fidelity the flow around 
the rotor. In our article we proposed a modeling method and a simulation of 
the forces applied to the horizontal axis wind rotors turbines using the 


Modeling application of the blade elements method to model the rotor and the vortex 
Vortex method method of free wake modeling in order to develop a rotor model, which can 
Wind power be used to study wind farms. This model is intended to speed up the 
calculation, guaranteeing a good representation of the aerodynamic loads 
exerted by the wind. 
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1. INTRODUCTION 

Use of wind energy as a source of renewable energy was considerably expanded in the 20th century, 
with the growing penetration of wind energy generation into power grids [1], [2]. Lanchester was the primary 
predicter of a perfect wind’s greatest power [3]. One of the faremost vital challenges within the aerodynamics 
of wind turbines is to predict precisely the forces on the blade, which are modeled by various approaches 
with varying degrees of complexity and precision, including the method of the blade element the vortex 
model and approaching the Navier-Stokes equations using CFD. Today, a BEM-based engineering model is 
commonly used to analyze the aerodynamic performance of a wind turbine where it relies on a normal and 
uniform flow hypothesis and on aerodynamic charges acting on an actuating drive rather than on a finite 
number of blades [4]. Hansen, et al., [5] show that the accuracy of the BEM under normal load conditions is 
reasonable if sufficient aerodynamic profile data is given. The method of the blade element is defined as the 
improved model of the Rankine-Froude dynamics theory [6], [7], and that was the firs model to predict the 
influx rotor speed where it is assumed that the rotor can be replaced with uniformly charged disk. CFD 
method is to resolve the differential equations of fluid mechanics. Navier-Stokes equations or simplification 
as the Euler equations, potential, potential linearized and RANS, by algorithms and numerical methods [8]. 
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Comprehensive review is presented in Sanderse, et al., [9] on the use of CFD codes for wind turbine 
calculations. The analytical techniques are easily defined according to the turbulence modeling method used, 
which is of particular significance for the modelling of wind turbine wakes. The vortex methods are another 
approach, which is more reliable than the theory of momentum, and less expensive than CFD. Methods based on 
vorticity have various formulations, from basic analytical models to more complex numerical methods [10]. The 
vortex theory can also be used to estimate the aerodynamic efficiency of wind turbines, based on the potential 
flow, non-viscous. The advantage of the method is that it reduces the complexity of the flow and it has 
generally been used for airfield and aircraft aerodynamic analysis. Despite not being able to predict the 
viscous phenomena such as separation and boundary layer, drag, its combination with aerodynamic profile 
data actually tabulated a powerful method for predicting the fluid flow. Methods of free vortex wake have 
developed as robust, flexible tools for modeling the aerodynamic loads [11]. Yu, et al., [12] developed a free 
wake method that uses a vortex ring model with a semi-infinite cylindrical vortex tube where near wake is 
simplified in to a series of vortex rings and the far wake with a semi-infinite cylindrical vortex tubes. 
Branlard and Gaunaa [13] developed approach based on the use of vortex-based methods including 
cylindrical and helical wake models to study the unstable effects of yaw, shear and turbulent flows. Another 
work of Vaal [14] presented an aerodynamic model that combined Kutta—Joukowski theory with the blade 
element model, which predict the blade charges accurately and effectively using the time-evolution of the 
induced speeds. Recent work by Jing Dong, et al., [15] suggested a modified method of free-wake vortex rings to 
resolve the aerodynamic load around wind turbines with a horizontal axis, on fixed and floating support 
structures. 

In this work, we propose an aerodynamic model which allows to quickly simulate the flow of a wind 
turbine with horizontal axis. This model makes it possible to precisely determine and calculate the position of 
the vortices and the induced radial and axial speeds at any point in the wake. A free wake vortex method has 
been implemented coupled with the blade elements method which allows us to quickly determine the 
aerodynamic loads. 


2. PRESENTATION OF THE MODEL 

During operation of a wind, the helical wake is formed behind a rotor, the latter which is modeled 
with a series of vortex rings [16]-[18]. For the modeling of the wind turbine rotor, one can use the theory of 
blade element which consists the blade divided into several parts. 


2.1. Modeling of the wake 

The wake is modeled by a series of vortex rings Figure 1, ready to move axially and to widen 
radially, able to interact with the impact of the speed field induced in the wake. Increased passage of the 
blade refers to the separation of the intensity of the vortex ring and the radius R, which is proportional to the 
axial force applied to the rotor. 





(b) 





Figure 1. (a) Wind wake visualization [19], (b) representation of the helical track by vortex rings 


1.2. Velocity induced by a vortex ring 
The velocity ild induced with the vortex rings is determineted using the law of Biot-Savart. Figure 2 
shows geometry of vortex rings with a radius R, located in plane xy. The induced velocity is calculated as (1): 
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I is the intensity of the vortex segment. 

Modeling the vortex rings by law Biot-Savart involves singularities within the vortex cores, so the 
obtained solution is, singular if the analysis point P falls on the ring, to avoid singular behavior, a 
regularization parameter 6 is introduced. Figure 3 depicts the axial velocity induced with and without the 
regularisation parameter. 


Induced axial velocity vz 


® p(x,y.zZ) 





Adimensional distance, r/R 


Figure 2. Geometry of the vortex ring in the plane xy Figure 3. Axial velocity induced with and without 
the regularisation parameter 
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Additionally, the precise resolution is voiced with the elliptical integrals of first and second kind at 


any point in the flux field [20] Thus, the radial and axial parts of the induced velocity are expressed in the 
following (3) (4) Vz Vr for the ring vortex with the radius R and also the intensity I. 
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Where the elliptic integrals K(s) and E(s) complete first and second kind in parameter S (5): 


fan (5) 
Z +(R-rf +5” 


The axial location of the vortex ring in the free wake is calculated with the addition of a wind speed U5, 
the self-induced velocity and mutually induced velocities of control point, as suggested by Heemst, et al., (6) [21]. 


j 
V 7 Uy T > v T V out (6) 


i=] 


Self-induced velocity is defined at the point situated in the vortex ring for radius R and radius 
core 7, according to the (7), as suggested by Lewis [22]. 
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Vaui = = (ve ee ~~ r) (7) 
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The radial velocity of the vortex ring is determined in a control point to the periphery of the vortex 
ring by adding all mutually induced velocities (8). 


j 
= DVM aut (8) 


i=l 


From Figure 4 We see that the vortex rings, which are formed right behind the rotor, are located in 
an area of high radial velocity, which causes their diameter to increase rapidly when away from the rotor. 
Such radial velocities, measured by the control point at the edge of the vortex ring, are reduced to zero. 
Figure 5 demonstrate the axial velocity as a function of the distance z of the rotor plane on the edge of the 
rings. It is remarkable that the vortex ring behind the rotor travels initially more slowly than those located 
further away from the rotor plane. 





z [m] z [m] 
Figure 4. Distribution of radial velocity of the vortex Figure 5. Distribution of axial velocity along the 
ring along the wake wake 


2.3. Modeling of the blades 

The aerodynamic modeling of the blade is of great importance in the design of the wind rotors, 
because very detailed calculations have to be carried out in order to determine the different essential 
dimensions of the blade. Use the blade element method to determine the force applied to the rotor blade, 
Figure 6 illustrate the speed triangle of blade element. In this approach, the blade is divided into segments of 
k, or each is assumed to be independent of the others. For each section the induced velocity is evaluated 
axially and tangentially by the vortex system, as shown by Karpatne, et al., [23]. 
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Figure 6. Speed triangle of blade element 
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The angle of incidence is given by (9) and (10): 


&=p-p (9) 
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The coefficient of local lift cl (1) is determinated for each ai after consulting the data table of the 
airfoil of the blade. The total thrust T generated by the rotor is given by (11): 


T= = 5 New’ cl(i)e(i)Ar (11) 


In order to calculate the intensity for each blade, the Kutta-Joukowski theorem is used. The 
equivalente vortex intensity can be determined from the total turbine thrust (12) [24] 


T= Tef Qrdr (12) 


By assuming that all the resulting flow is dragged from the tip of the blades and adding the total 
thrust of all the blades, the intensity of the vortex ring is determinated by (13) [24]. 


2T 
a (13) 
PNR*Q 
Where P is the circulation that will be released from the vortex ring. This iterative method for I’ is performed 
at each time step At, as shown in Karpatne, et al., (14) [23]. 


2T 


At = — 
NQ 


(14) 


2.4. Scheme of temporal integration 
At each time step, the total field of velocity induced by all the rings shall be calculated. For the 

passing of time a predictor-corrector scheme has been selected as shown in (15) and (16) [25]. As shown in 

(15) and (16) are used for the first iteration to predict a solution of axial and radial positions of the vortex 

rings. 

r =r" + Atv, (7,” P zr) 
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The Adams-Bashforth method is used for precision and stability reasons to determine the location of 
the vortex rings. As shown in (17) to (18). 


r 1 
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Thereafter, the speeds ((vi-(rř, Z; ), Viz(1%j Z; )) are determinated in this new position). The ring 
location is determined by average intermediate velocity and new velocies (19) and (20). 
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3. RESULTS AND DISCUSSION 

In this work, we have used two types of wind turbines to study and simulate the forces applied to 
wind turbine rotors. The first case is a wind-rotor type NRELSMW at wind speed of 8 m/s. The NREL report [26] 
contains some basic parameters for turbines, aerodynamic profiles, aerodynamic characteristics and other 
turbine details used in this model. Figure 7 shows the distribution of vortex rings in the wake, the rings move 
axially and expand radially. 

This numerical simulation us allows the following results presented in Figures 9 and 10 which show 
the distribution of the intensity I’, the tangential and the axial forces distribution along the blade. In Figure 8 
it is noted that the intensity I blade tip is 39 m7/s this intensity increases proportionally to the value of 0.87 
r/R and decreases rapidly in the blade root. 

The tangential and the axial forces distribution in the rotor plane is sown in Figure 9 (a) (b). Note 
that these forces increases slightly at the end of the blade and stabilizes, after these forces decreases rapidly at 
the blade roots. The second case is a wind turbine: Gamesa G58 with an NACA 4412 profile 3 blades, the 
rotor at a radius of 58 meters from the results obtained in Figure 10 shows the power recovered by the turbine 
Gamesa according to the wind speed. 
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Figure 7. Vortex rings distrbution in the wake Figure 8. Distribution of intensity I’ along the blade 
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Figure 9. (a) The tangential forces distribution in the rotor plane, (b) the axial forces distribution in the rotor 
plane 
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Figure 10. Power produced by gamesa G58 according to the wind speedy 


CONCLUSION 
In this paper, we presented a work that involves the modeling study of the forces applied to the wind 


turbine rotor axis horizontal. To do this, we have proposed a model that combines the method of blade 
element for determination of aerodynamic forces and vortex theory to calculate the velocity field induced in 
the rotor plane. This model has allowed the implementation of a model digital capable of performing 
simulations with acceptable accuracy and a time adequate calculation. Numerical simulations are performed 
and the simulation results of the two cases of wind turbines show the validation of this model to determine 
the wind turbine performance and various aerodynamic parameters blades. It is possible to consider the 
proposed model as operational and can also be used to simulate unstable cases. 
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